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ABSTRACT 

Observations of the Vela pulsar region with the Chandra X-ray observatory 
have revealed the fine structure of its synchrotron pulsar- wind nebula (PWN), 
which showed an overall similarity with the Crab PWN. However, contrary to the 
Crab, no firm detection of the Vela PWN in optical has been reported yet. To 
search for the optical counterpart of the X-ray PWN, we analyzed deep optical 
observations performed with different telescopes. We compared the optical im- 
ages with those obtained with the Chandra ACIS to search for extended emission 
patterns which could be identified as counterparts of the X-ray nebula elements. 
Although some features are seen in the optical images, we find no correlation with 
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the X-ray structure. Thus, we conclude that the diffuse optical emission is more 
likely associated with filaments in the host Vela SNR. The derived upper limits 
on the optical flux from the PWN are compatibile, within the uncertainties, with 
the values expected on the basis of the extrapolations of the X-ray data. 

1. Introduction 

A compact (~ 4') X-ray nebula around the Vela pulsar was first detected in soft X-rays 
by Harnden et al. (1985) with the Einstein High Resolution Imager (HRI). These authors 
suggested that it is a pulsar-wind nebula (PWN) powered by relativistic particles ejected by 
the pulsar. The soft X-ray emission from the Vela X-ray nebula was further studied with 
ROSAT (Ogelman, Finley & Zimmerman 1993; Markwardt & Ogelman 1998, and references 
therein). Based on its "kidney-bean" shape, Markwardt & Ogelman (1998) proposed an 
interpretation in terms of a bow-shock produced by the supersonic motion of the pulsar 
through the ambient medium. 

The Vela pulsar and its PWN have been later observed with both the High Resolu- 
tion Camera (HRC) and the Advanced CCD Imaging Spectrometer (ACIS) on board the 
Chandra X-ray observatory (Pavlov et al. 2000, 2001a,b, 2003; Helfand, Gotthelf, & Halpern 
2001; Kargaltsev et al. 2002). Thanks to the excellent angular resolution of Chandra, the 
morphology of the PWN was resolved in a complex structure resembling that of the Crab 
PWN. Such a structure cannot be explained by a simple bow-shock model. The brighter 
part of the PWN, with a size of about 2', shows an approximately axisymmetric structure, 
with two arcs, a jet, and a counter-jet, embedded into an extended diffuse emission. The axis 
of symmetry, which can be associated with the pulsar rotational axis (Pavlov et al. 2000), 
coincides, within the errors, with the direction of the pulsar's proper motion (P.A. = 301°— 
e.g., Caraveo et al. 2001a). 

The overall spectrum of the PWN can be described by a power law with an average 
spectral (energy) index a ~ 0.5 (see Kargaltsev et al. 2002 for details), which can be inter- 
preted as synchrotron emission of relativistic electrons and/or positrons. The same electron 
distributions should emit optical and radio synchrotron radiation, provided that the electron 
power-law spectrum extends down to sufficiently low energies, < 300isl{ 2 BZ^ 2 GeV (where 
v = 10 14 ^i4 Hz is the radiation frequency, and B = 10~ 4 5_ 4 G is the magnetic field), and 



1 Based on observations collected at the European Southern Observatory, La Silla, Chile 

2 Bascd on observations with the NASA/ESA Hubble Space Telescope, obtained at the Space Telescope 
Science Institute, which is operated by AURA, Inc. under contract No NAS 5-26555 
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synchrotron self-absorption plays no role. 

Indeed, highly polarized 60% at 5.2 GHz) extended radio emission has been recently 
detected around the Vela pulsar (Lewis et al. 2002; Dodson et al. 2003), covering a region ~4 
times larger than the X-ray PWN, as observed in the case of the Crab. Most of the diffuse 
radio emission comes from two lobes (see Figs. 2-4 in Dodson et al. 2003) — the southeast 
lobe of an area of 18 arcmin 2 (F„ = 760 ± 100 mJy at 5.2 GHz) and a brighter northeast lobe 
of an area of 5.3 arcmin 2 (F u = 290 ±50 mJy at 5.2 GHz). Search for radio emission from the 
X-ray-bright compact nebula was hampered by the brightness of the pulsar. Although some 
emission was detected (e.g, about 30 mJy at v — 2A GHz, in a 0.73 arcmin 2 area around 
the pulsar), it can well be an unsubtracted pulsar contribution and should be considered as 
an upper limit on the compact PWN emission. 

One could expect to see optical nebula around the Vela pulsar. Indeed, optical PWNe 
have been observed in several other young pulsars. The most famous example is the Crab 
PWN (see, e.g., Hester et al. 2002, and references therein) which exhibits a complex structure 
featuring bright knots and sharp arc-like features dubbed "wisps" (Scargle 1969) located at 
the inner boundary of the X-ray torus and probably associated with a termination shock in 
the pulsar wind (e.g., Hester et al. 1995; Kennel & Coroniti 1984). Also, the comparison 
between the Hubble Space Telescope (HST) and Chandra observations of the LMC pulsar 
B0540-69 has clearly shown the presence of an optical PWN which follows exactly the 
morphology of the X-ray PWN (Caraveo et al. 2001b, 2003). 

Contrary to the Crab and PSR B0540-69, searches for the Vela PWN in optical have 
been inconclusive so far. The only marginal detection of an optical PWN was reported by 
Ogelman et al. (1989), who observed the Vela pulsar field in the V and B bands with the 
ESO 2.2 m telescope. These authors found evidence for optical diffuse emission (typical size 
~ 2'), with an average surface brightness of about 26 mag arcsec -2 . However, the presence 
of bright filaments from the host SNR, as well as of several bright stars in the field, made it 
difficult to assess the reality of the putative optical nebula. 

The increase in sensitivity and angular resolution provided by the HST prompted us 
to carry out a new optical investigation of the Vela pulsar region. In the following, we 
discuss spatial correlations between Chandra ACIS images (Pavlov et al. 2001b; Pavlov et 
al. 2003) and the optical ones collected by both the HST WFPC2 (Mignani and Caraveo 
2001; Caraveo et al. 2001a) and ESO NTT (Nasuti et al. 1997) and VLT 
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2. X-ray and optical observations 

2.1. Chandra observations 

The multiple observations of the Vela PWN with the Chandra ACIS are described 
by Pavlov et al. (2003). In all the observations the PWN was imaged on the same back- 
illuminated chip S3, more sensitive to soft X-rays. Although the individual images show 
some variability of the PWN elements, the overall PWN morphology remains stable. To 
obtain a very deep X-ray image to be compared with the optical ones, we stacked seven 
consecutive individual images, collected from 2001 November 25 through 2002 April 3 for a 
total exposure time of 141 ks, and adaptively smoothed the combined image. The accuracy 
of the image co-alignment was about &!7. The smoothed X-ray image of the PWN, in the 
energy band 1-8 keV, is shown in Figure 1. The logarithmic brightness scale was chosen 
to increase the dynamic range. The jet and the counter-jet, northwest and southeast of 
the pulsar, respectively, and two arc-like structures (inner and outer arcs) are the dominant 
features of the bright inner nebula (upper panel). The inner nebula is surrounded by a 
bean-shaped region of diffuse emission (outer nebula), ~ 2' x 2' in size (lower panel). An 
elongated region of fainter extended emission southwest of the bright inner PWN is seen up 
to the edge of the ACIS image (~ 4' from the pulsar). A 100"-long outer jet is seen northwest 
of the inner PWN. The spectrum of the bright inner PWN, within a 45" radius around the 
pulsar, can be fitted by a single power-law model with a spectral index a = 0.47 ± 0.02 and 
unabsorbed flux of 9.1 ±0.7 x 10~ n ergs s _1 cm" 2 , in the 0.1-10 keV range (Kargaltsev et al. 
2002). The unabsorbed surface brightnesses for the outer and inner arcs are about 6 x 10~ 14 
and 3 x 10~ 14 ergs s _1 cm~ 2 arcsec~ 2 , while the spectral indices a are about 0.4 and 0.3, 
respectively. 

2.2. Optical wide-band imaging 

Optical observations of the Vela pulsar field were collected at different epochs using 
different telescopes, instrumentation, observational set-ups, and filters. In addition to the 
recent high-resolution imaging data obtained with the HST, our data set includes archived 
ground-based images obtained from ESO with the NTT and the VLT. Table 1 gives a detailed 
summary of the available observations. 
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2.2.1. HST 

HST observations of the Vela pulsar field were collected between 1997 June and 2000 
March with the WFPC2. Several images have been taken through the filters 675W (A = 
6717A; AA = 1536A), 8UW (A = 7995A; AA = 1292A) and 555VF (A = 5500A; AA = 
1200A) with exposure times ranging between 2000 s and 2600 s. These data were origi- 
nally taken as part of two independent programs aimed at the study of the multicolor flux 
distribution of the pulsar (Mignani & Caraveo 2001) and at the measure of its parallactic 
displacement (Caraveo et al. 2001a). To achieve the maximum spatial resolution, in all the 
exposures the pulsar was located near the center of the Planetary Camera chip, with a pixel 
size of 0"045 and a field of view of 35" x 35". The standard HST pipeline was applied for 
the image reduction and for the photometric calibrations. Co-aligned exposures were finally 
combined to filter out cosmic ray hits (see Mignani & Caraveo 2001 and Caraveo et al. 2001a 
for a detailed description of the observations and data reduction). 

For each observation, the IRAF/STSDAS task wmosaic was used to obtain a mosaic 
of the four WFPC2 images and correct for the instrumental geometric distortion of the 
cameras (Casertano & Wiggs 2001). To increase the signal-to-noise ratio, we combined all 
the available 555VF filter images. Since the observations at different epochs were performed 
with different telescope pointings and roll angles, to superpose the frames we had to follow 
the procedure successfully applied in previous astrometric works, using a common reference 
frame defined by a grid of reference stars' coordinates (see Caraveo et al. 2001a and references 
therein for further details). Images were registered within 0.1 WFC pixels (~ 10 mas), an 
accuracy adequate to achieve the goals of the present investigation, co-aligned and stacked. 

2.2.2. NTT 

The Vela pulsar field was observed with the NTT at the La Silla observatory in January 
1995 (Nasuti et al. 1997). To achieve the highest sensitivity both at longer and shorter 
wavelengths, the EMMI camera was used in both its Red (EMMI-R) and Blue (EMMI-B) 
arm modes. The camera is a TEK 2048 x 2048 pixels CCD with a pixel size of 0"27 (field of 
view 9!2 x 8'.Q ) and 0"37 (field of view Q'.2 x Q'.2) in the EMMI-R and EMMI-B configuration 
modes, respectively. Images were collected in the Johnson's U (A = 3542A; AA = 542A), B 
(A = 4223A; AA = 941 A), V (A = 5426A; AA = 1044A), and R (A = 6410A; AA = 1540A) 
bands with exposure times between 15 and 40 minutes. Data reduction and photometric 
calibration were carried out as described in Nasuti et al. (1997). 
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2.2.3. VLT 

Observations of the field were performed in April 1999 with the first Unit Telescope 
(Antu) of the ESO VLT located at the Paranal Observatory. Images were obtained using 
the FOcal Reducer and Spectrograph #1 (FORS1) instrument, a four-port 2048 x 2048 
CCD detector which can be used both as a high/low resolution spectrograph and an imaging 
camera. The instrument was operated in imaging mode with a 1 x 1 binning and at its 
standard angular resolution of 0.2 arcsec/pixel, with a corresponding field of view of 6(8 x 6(8. 
Two 300 s exposures were taken in the Bessel filters R (A = 6570A; AA = 1500A) and 
I (A = 7680A; AA = 1380A) under good seeing conditions (< 1"). The images have 
been retrieved from the ESO public archive and reduced following standard procedures for 
debiassing and flatfielding. Since no standard stars were observed, absolute flux calibration 
was performed using as a reference a set of secondary stars detected in the same passbands 
in the NTT images. 



2.3. Narrow-band imaging 

Three 20-min exposures of the Vela pulsar field were taken in H a (A = 6588.27A, 
AA = 74.31A) on 1999 April 4 with the Wide Field Imager (WFI) at the ESO/MPG 2.2 m 
telescope. The WFI is a wide field mosaic camera, composed of eight 2048 x 4096 pixel 
CCDs, with a scale of 0.238 arcsec/pixel, providing a field of view of 33(7 x 32(7. To 
compensate for the loss of signal due to the interchip gaps (23'.'8 and 14'.'3 along right ascension 
and declination, respectively), the second and third exposures were taken with a relative 
offset of 30" in RA and Dec. The images were taken in fairly good seeing conditions with 
FWHM~ 0'.'9. Data reduction with the usual CCD processing steps of bias subtraction, 
flat-fielding and trimming of images, was performed with IRAF 3 using the Mosaic CCD 
reduction package (MSCRED). Individual exposures have been finally coadded and cleaned 
of cosmic-ray hits using the DRIZZLE software (Hook & Fruchter 1997) and the final image 
has been corrected for the effects of fringing that affects the CCDs in the red part of the 
spectrum (> 6500 A). 



3 IRAF is distributed by NO AO, which is operated by the AURA, Inc., under cooperative agreement with 
the NSF. 
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2.4. Image superposition: X-ray versus optical 

The direct comparison of images taken in different energy bands (e.g., optical and X- 
rays) can be achieved through accurate image superposition. Following the approach used 
by Caraveo et al. (200fb), we superimposed the X-ray Chandra image onto the optical ones 
with respect to the absolute (a, 6) reference frame, relying on the astrometric solution of 
each image. We used the combined Chandra ACIS image, where individual images were 
co-aligned with an accuracy of Of! 7 (Pavlov et al. 2003). The error of the absolute Chandra 
aspect solution 4 is typically about Of! 6. Therefore, we adopt 1" as a reasonable estimate for 
the uncertainty of the absolute X-ray astrometry. 

For the WFPC2 images, the default astrometric solution across the focal plane is derived 
from the coordinates of the two guide stars used to point the telescope, which are taken as 
a reference to compute the astrometric reference point and the telescope roll angle. The 
accuracy of the HST astrometry is limited by the intrinsic error on the absolute coordinates 
of the GSC1.1 (Guide Star Catalog) stars (Lasker et al. 1990) which are used for the telescope 
pointing. According to the current estimates, the mean uncertainty of the absolute positions 
quoted in the GSC1.1 is about 0'.'8 per coordinate (Biretta et al. 2002). For this reason, we 
have recomputed a new astrometric solution for the WFPC2 images by using as a reference 
the positions of several stars from the USNO-A2.0 catalog (Monet 1998) identified in the 
field of view. The ASTROM software 5 was used to compute the pixel-to-sky coordinates 
transformation. The rms of the astrometric fit was ~ / .'08, per coordinate. 

The astrometric calibrations of the the ground-based optical images was computed as 
described above. The astrometric fits yielded rms values of ~ 0'.'13 and ~ 0'.'22 per coordi- 
nate for the VLT and NTT images, respectively. For the WFI image, the astrometry was 
computed separately for each of the eight CCD chips, yielding an average rms of ~ 0'.'3. 

In all cases, the final uncertainty on our astrometric solutions was evaluated by ac- 
counting both for the rms of the astrometric fits and for the propagation of the error on the 
intrinsic absolute accuracy of the USNO-A2.0 coordinates, which is of the order of / . / 2-0'. / 3. 

The absolute frame registration between the optical and X-ray images turned out to 
be accurate within ps 1", i.e., compatible with the overall uncertainties of the absolute 
astrometry of each frame. Unfortunately, the strong pile-up of the Vela pulsar in the ACIS 
images hampered the use of the pulsar position for better co-alignment between the X-ray 



4 see http : // cxc . harvard. edu/cal/ASPECT/celmon/ 

5 http://star-www.rl.ac.uk/ Software / soft ware . ht m 
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and optical images. 

3. Results 

3.1. Search for a compact optical nebula 

As a first step, the central part of the X-ray PWN field (corresponding to the inner 
nebula and the outer nebula) was inspected to search for optical counterparts of the complex 
structures seen with Chandra. Our starting point was the combined WFPC2 image, 
which is by far the deepest optical image of the Vela pulsar field. The final image is shown 
in Figure 2, where we superimposed the X-ray contour map obtained from the combined 
Chandra ACIS exposure of the region. Although a number of complicated patterns of diffuse 
emission are present in the ~ 3' x 3' field of view, no optical counterparts of the X-ray features 
seen in Figure 1 can be identified, nor any other structure symmetric with respect to the 
axis of symmetry of the X-ray PWN. We thus conclude that the PWN is undetected in the 
optical. 

To compute the upper limit on the optical surface brightness of the PWN, an accurate 
mapping of the background is required. This is complicated by the presence of diffuse, non- 
uniform, emission patterns which show sharp surface brightness variations on angular scales 
as small as 5 arcsec. For this reason, we evaluated the background level in 200 cells of 
1 arcsec 2 each, selected in a number of star-free regions across the whole image. Statistical 
errors on the number of counts per cell were of the order of 0.3%. The background level, 
~ 22.3 ST magnitudes 6 arcsec" 2 on average, was found to vary typically by 4%-5% across 
the whole field, with a maximum variation of ~ 7%. 

For an extended source, the upper limits on the surface brightness scale with the de- 
tection area A as A^ 1 ^ 2 . An optimally chosen area should be large enough to reduce the 
statistical errors due to the background fluctuations but it should be smaller than the scale 
of sharp background variations. According to our mapping of the background, a detection 
area of 10 arcsec 2 represents a resonable optimization. 

For the area chosen, the measured background variations affect the upper limit on the 
flux of an extended source (at a 3 a level) by no more than p^3%, corresponding to surface 
brightness variations below 0.1 magnitudes arcsec" 2 . 

While variations in the sky background across the image play a minor role, we found 



6 http://www. stsci.edu/instruments/wfpc2/wfpc2_doc. html 
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that the derived surface brightness upper limits are different in different regions of the image. 
This is mainly due to the non-uniform coverage of the field performed by the WFPC2. First, 
because of the intrinsic differences in the pixel size and in the physical characteristics of the 
CCDs, the PC and the WFC chips contribute differently to the instrumental background 
per unit area. In addition, since the five WFPC2 observations listed in Table 1 were 

executed at different epochs and with different telescope roll angles, the exposure map varies 
across the field. For instance, the central region of the field, which is covered by the PC 
(35" x 35"), reaches an integration time of 11800 s, while the outer regions, covered by the 
three WFC chips, have integration times varying between 4 600 and 11800 s. 

Both effects clearly affect the evaluation of the surface brightness upper limits in different 
regions of the PWN, as they are covered differently by the four WFPC2 chips. As it is seen 
from Figure 2, some regions of the inner nebula (inner arc, jet, and counter-jet) fall entirely 
within the PC, while the outer arc is coincident with the inter-chip gaps and is covered 
partially by the PC and partially by the WFC chips. On the other hand, the outer nebula 
is entirely covered by the WFC chips. 

Taking all these effects into account, we have computed the 3a upper limits on the 
optical surface brightness of both the inner nebula (inner arc, jet, counter-jet, outer arc) 
and the outer nebula. We note that although our upper limits have been computed for a 
detection area of 10 arcsec 2 , they can be easily rescaled to any other area. 

Using the HST pipeline photometric calibration, we computed an upper limit of 28.1 and 
28.0-28.5 ST magnitudes arcsec -2 for the inner nebula and for the outer nebula, respectively. 
As we mentioned before, because of the non-uniformity of the exposure map the upper limit 
for the outer nebula turned out to be slightly position-dependent. To correct these limits 
for the interstellar reddening, we use the extinction Ay = 0.2, consistent with the hydrogen 
column density estimated from the X-ray observations of the Vela pulsar (Pavlov et al. 
2001a). The corrected limits are 27.9 and 27.8-28.3 ST magnitudes arcsec -2 . 

The same analysis was then repeated for the other available WFPC2 images (675W and 
814VF) as well as for the NTT (UBVK) and VLT (VI) ones (see Table 1) but no evidence 
for a compact optical nebula was found in either of these datasets. The computed upper 
limits are summarized in Table 2. 



3.2. Search for an extended nebula 

Since the emission from the PWN could, in principle, be visible at optical wavelenghts 
on larger angular scales with respect to the X-rays, as it has been observed in radio (see 
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§1), we took advantage of the larger field of view provided by the NTT images to search 
for extended features up to distances of 3 arcmin. In particular, we searched for diffuse 
optical emission at the position of the southwest extension of the X-ray nebula (see Fig. 1, 
lower panel). To go as deep as possible, we have combined all the available NTT UBVR 
band images (Fig. 3, upper panel). The combined image shows many different enhancements 
in the background, with a rather complex spatial distribution. However, none of them can 
be firmly correlated with the known X-ray features. In addition, we note that almost all 
the diffuse emission patterns seen in the NTT UBVR image can be also identified in the 
ESO/2.2m H a (Fig. 3, lower panel). This suggests that they are most likely associated with 
the bright filaments of the Vela supernova remnant. We note that the computed NTT upper 
limits on the optical emission of the compact X-ray nebula can be applied also at larger 
distances from the pulsar. The maximum variation in surface brightness for source-free 
regions, due to the complicated distribution of diffuse emission, is found to be of order 7% 
even in the outer regions of the NTT field, which corresponds to < 0.1 magnitudes arcsec -2 
variations of the upper limit across the field (see §3.1). We note that for the NTT images 
the measured upper limits in the U and B bands do not apply to the region ps 2.5 arcmin 
southwest of the pulsar position, close to the southwest edge of the extended X-ray nebula. 
The background in this region is unrecoverably polluted by the presence of a bright (B rs9) 
B star that is not visible in the V and R band images because of the slightly narrower field 
of view. 



4. Discussion 

It is interesting to compare the measured upper limits on the optical brightness of the 
Vela PWN with the X-ray and radio observations. The observed X-ray spectra of the Vela 
PWN are described by a power law with spectral indices a ~ 0.3-0.4 and a ~ 0.4-0.5, for 
the arcs and the extended diffuse emission southwest of the pulsar, respectively (Kargaltsev 
et al. 2002). The X-ray and radio surface brightness spectra, together with the optical 
upper limits, are shown in Figure 4 for the outer arc (upper panel) and the diffuse emission 
southwest of the pulsar (lower panel). Since the inner PWN structures were not resolved in 
the radio (see §1), only radio upper limits are plotted in the upper panel. From the plot we 
see that both the optical and radio upper limits for the outer arc are within the uncertainty 
of the extrapolation of the X-ray spectrum towards lower frequencies. 

In the case of the southwest diffuse emission, we note that the X-ray-to-radio extrapola- 
tion is well below the measured radio brightness values. This apparent inconsistency can be 
explained by the fact that the radio and X-ray brightnesses were measured in different areas 
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(the X-ray image is substantially smaller than the radio image). Moreover, even within the 
smaller X-ray field-of-view, it is seen that the X-ray brightness of the southwest diffuse emis- 
sion fades towards the region of maximum radio brightness (which is outside the X-ray field 
of view). Such behavior can be explained by radiative (synchrotron) and adiabatic cooling 
of the expanding cloud of relativistic electrons, which can result not only in the increase 
of the spectral index, but also in the shift of the lower-energy boundary of the power-law 
spectrum towards lower energies. Therefore, we can expect maximum optical brightness to 
be observed in between the regions of the X-ray and radio maximum brightnesses. Although 
the X-ray and radio spectra in the lower panel of Figure 4 cannot be directly compared with 
each other, a crude estimate of the expected optical brightness can be obtained by connect- 
ing the radio and X-ray points. This yields ~ 3 x 10~ 32 erg cm -2 s _1 Hz -1 arcsec -2 , i.e., 
about 3 magnitudes below our measured optical upper limit. 

To summarize, we conclude that it is very likely that just a slightly deeper optical 
observation would allow one to detect the outer/inner arc (and other bright elements of the 
inner PWN), while much deeper optical observations are needed to detect the emission from 
the Vela PWN at large. However, owing to the presence of background/foreground emission 
from the SNR as well as numerous bright stars in the field, longer exposures would not help to 
detect the optical PWN — neither its X-ray-bright central part nor the radio-bright outskirts. 
To get rid of the contaminating emission, one should observe the field at UV wavelengths, 
where the radiation of most of the field stars would be much dimmer. Observing in polarized 
light appears even more promising. Since the synchrotron emission from the PWN should be 
highly polarized (as confirmed by the radio observations), polarimetry observations should 
allow one to minimize the contamination from field sources and provide a clean PWN image. 
Polarimetry observations of the Vela pulsar field have been recently obtained by Wagner 
and Seifert (2000) but they did not yield a conclusive result. Deeper, higher-resolution, 
observations are needed to unveal extended polarized emission from the PWN. 
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Fig. 1. — : X-ray images of the Vela PWN, combined of seven individual Chandra ACIS images (see 
Pavlov et al. 2003). The upper panel (l'.QO x lf35) shows the bright inner PWN around the Vela 
pulsar (marked by the cross), with the outer and inner arcs, northwest inner jet, and southeast 
counter-jet. The brightness scale in the lower panel (6'.3 x 5'A) was chosen so as to show the faint 
nebular structures — e.g., the northwest outer jet and the dim extended emission southwest of the 
bright inner PWN. North is up, East to the left. 



Fig. 2. — : Image of the Vela pulsar field obtained by combining all the WFPC2 observations 
listed in Table 1 (North to the top, East to the left). The gaps among different CCD chips are 
evident. The overall integration time on the central part, corresponding to the PC field of view, is 
11800 s (see text). The overlayed contours (logarithmic scale) correspond to the X-ray intensity 
maps obtained from the Chandra ACIS image of the field in the energy band 1-8 keV. The point 
source within the innermost X-ray contour is the optical counterpart of the Vela pulsar. 
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Fig. 3. — : The upper panel shows the combined UBVR image of the Vela pulsar field obtained 
from the JVTT/EMMI observations listed in Table 1. The lower panel shows the ESO/2.2m H a 
image. In both cases the image size is ~ 4' x 4' . North is up, East to the left. The X-ray contour 
plots of the Vela PWN derived from the seven ACIS observations (§2.1) are overlayed. 
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Fig. 4. — : Spectra of surface brightness (in erg cm -2 s _1 Hz -1 arcsec -2 ) in X-rays (solid lines) and 
radio (points), together with the optical upper limits, for the inner/outer arc (upper panel) and 
diffuse emission southwest of the pulsar (lower panel). Expected brightness levels in optical, based 
on extrapolations of the X-ray and radio data, are shown with dashed lines. 
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Table 1. Available optical datasets for the Vela pulsar field. 



Date 


Telescope 


Instr. 


Filter 


A (AA) 


Exp. 


Ref. 


Tan 1 QQS 


NTT 

1 \ J. J. 


FMMT-R 

J_J1V J. IV J.X 1 J 


TJ 


-.7 '-J T: Z. i X. 1 UliiA J 


4800 


(1) 


Jan 1995 


NTT 


EMM1-B 


B 


4223A (94lA) 


1800 


(1) 


Jan 1995 


NTT 


EMMI-R 


V 


5426A (1044A) 


1200 


(1) 


Jan 1995 


NTT 


EMMI-R 


R 


6410A (1540A) 


900 


(1) 


June 1997 


HST 


WFPC2 


555Vy 


5500A (1200A) 


2600 


(2) 


Jan 1998 


HST 


WFPC2 


555W 




2000 


(2) 


June 1999 


HST 


WFPC2 


555W 




2000 


(2) 


Jan 2000 


HST 


WFPC2 


hhhW 




2600 


(2) 


Jul 2000 


HST 


WFPC2 


555Vy 




2600 


(2) 


Mar 2000 


HST 


WFPC2 


6751^ 


6717A (1536A) 


2600 


(3) 


Mar 2000 


HST 


WFPC2 


8UW 


7995A (1292A) 


2600 


(3) 


Apr 1999 


VLT 


FORS1 


R 


6750A (1500A) 


300 




Apr 1999 


VLT 


FORS1 


I 


7680A (1380A) 


300 




Apr 1999 


2.2m 


WFI 


H a 


6588A (74.3A) 


3600 





Note. - First column lists the epoch of observation. Second and third 
columns show the telescope and the detector used for the observations, 
respectively. The filter names are listed in column four, with their pivot 
wavelengths and widths in column five. The total integration time per 
observation (in seconds) is given in column six. The last column provides 
the references: (1) Nasuti et al. 1997; (2) Caraveo et al. (2001a); (3) Mignani 
& Caraveo (2001). 
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Table 2. 3a upper limits to the surface optical brightness of the X-ray PWN structures. 



Observed Extinction-corrected 1 



Telescope 


Instrument 


Filter 


mag arcsec 2 


Flux 2 


mag arcsec 2 


Flux 2 


HST 


WFPC2 


555J¥ 


28.1 


0.21 


27.9 


0.25 








28.5-28.0 


0.15-0.23 


28.3-27.8 


0.18-0.28 




WFPC2 


Q75W 


27.5 


0.55 


27.3 


0.66 








27.9 


0.38 


27.7 


0.46 




WFPC2 


8UW 


27.7 


0.65 


27.6 


0.71 








28.1 


0.45 


28.0 


0.49 


NTT 


EMMI-B 


U 


26.4 


0.52 


26.1 


0.69 






B 


27.4 


0.47 


27.1 


0.62 




EMMI-R 


V 


27.1 


0.53 


26.9 


0.64 






R 


26.7 


0.59 


26.5 


0.71 


VLT 


F0RS1 


R 


27.0 


0.45 


26.8 


0.54 






I 


26.1 


0.81 


26.0 


0.90 



Hot A v = 0.2 

2 flux values are in units of 10 -30 ergs cm" 2 s _1 Hz -1 arcsec -2 

Note. - - The upper limits are computed for an area of 10 arcsec 2 . For the HST results, 
the first and second rows, for a given filter, are the upper limits measured in the PC and 
WFC chips, respectively. Due to the uneven exposure map of the combined WFPC2 555VT^ 
image (Fig. 2), slightly different upper limits are derived across the WFC field. For both 
the NTT and VLT, the upper limits apply to the overall X-ray nebula. 
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